A validated LC–MS/MS method for simultaneous determination of key glucocorticoids in animal hair for applications in conservation biology

A new method for the determination of main glucocorticoids (cortisol, cortisone, and corticosterone) in hair by liquid chromatography-tandem mass spectrometry was developed. Glucocorticoids were extracted from hair shafts using methanol followed by solid-phase extraction. A validation test was performed using hair from three species of wild mammals with different body size (0.2–800 kg), lifestyle (terrestrial, burrowing and arboreal species), social organization (living in herds or solitary), and different predicted type of hair glucocorticoids: European bison (Bison bonasus), European hamster (Cricetus cricetus), and Eurasian red squirrel (Sciurus vulgaris). Regardless of the species evaluated, the method shows good linearity for all analytes accompanied by satisfactory accuracy (91–114%) and precision (RSD < 13%). Depending on the analyte and hair origin, the calculated limits of quantification were between 0.05 and 1.19 ng/mL, which corresponds to 1.28–31.51 pg/mg. Using cortisol and cortisone as examples, we have demonstrated that measuring multiple glucocorticoids simultaneously provides more comprehensive information than solely concentrating on one, thereby contributing to a more balanced and reliable interpretation of the acquired results. However, the utility of cortisol metabolites as markers of stress response in keratinized tissues should be substantiated by additional experimental studies on targeted animals. We posit that this paper could serve as a crucial catalyst to prompt such experiments.

Monitoring of changes in levels of glucocorticoids in blood, saliva, or urine is a well-established strategy for linking psychological and health-related processes in medical, veterinary, and conservation studies [5][6][7] .Analysis of this kind of sample, however, allows for investigation of a time-point or short-term patterns of steroid hormones secretion and may be influenced by many concomitant factors, including the effect of sampling [8][9][10] .For this reason, matrices are also sought to determine the long-term, average glucocorticoids level of animals 11 .
Numerous studies have proven that keratinized tissues such as hair or claws can be useful matrices for the determination of glucocorticoids from the last weeks or months 9,[12][13][14] .It is known that glucocorticoids can diffuse from the bloodstream into the hair and then diffuse throughout the shaft but the kinetic of this process depends on the species basis and the thickness of the hair 15 .This variability makes it impossible to establish uniform experimental conditions for assessing the concentration of steroids in keratinized tissues across different animal species, and experimental protocols validated for one species may not yield reliable results for another.It should also be noted that the measurement of glucocorticoids can serve as a basis for implementing specific practical measures to improve the well-being of the studied population.Such circumstances arise, for example, in conservation biology, where the physiological state and welfare of animals form the foundation for implementing specific (often costly) conservation measures 16 .In this context, a reliable quantitative measurement of the long-term concentration of steroids in a given matrix is crucial.
Enzyme-linked immunosorbent assay (EIA) is frequently used to measure endogenous steroid contents in various biological samples due to its rapidness and simplicity 17,18 .However, this technique allows for the monitoring of a single analyte (e.g.cortisol), suffers from poor specificity due to cross-reactivity with other endogenous steroids and lipids and limited dynamic range 8,[19][20][21] .Paradoxically, the cross-reactivity of antibodies in EIA assays can be considered a useful attribute, since it enables the simultaneous detection of glucocorticoid metabolites with a similar chemical structure, which might be present in higher concentrations than the native (unmetabolized) compound within the sample.This approach, focusing on group determination of metabolites rather than native compounds, has found widespread application in studies examining fecal samples from various species of wild animals 22 .It is worth noting that selecting an inappropriate EIA for analysis may result in inaccurate representations of adrenocortical activity, as was for example observed in a study with hair samples 23 .
Besides EIA, numerous LC-MS/MS methods have been developed to measure steroid hormones (including glucocorticoids) in body fluids and hair extracts with good accuracy, precision, and selectivity 12,18,20,21,[24][25][26] .The LC-MS/MS measurement usually requires a small sample volume or amount and allows for the determination of numerous analytes in a single analytical run.However, LC-MS/MS analysis of keratinized tissue samples is known to be challenging due to low endogenous levels of steroid hormones and strong signal suppression caused by co-extracted interfering compounds (so-called matrix effect, ME).Thus, existing LC-MS/MS protocols were devoted to optimizing the extraction of steroid hormones from human 24,25 or animal hair samples 10 , but there are still significant methodological gaps that make research on different groups of animals difficult.These gaps relate in particular to the possibility of applying the achievements of modern analytics (commonly used, for example, in medicine) in the study of wild and endangered animal species in order to better understand their physiological condition and to make their conservation more effective 27 .
The overarching goal of this work was to develop and validate the new ultra-high performance liquid chromatography-electrospray ionization-tandem mass spectrometry (UHPLC-ESI-MS/MS) method for simultaneous quantification of cortisol, cortisone, and corticosterone in hair collected from wild animal species (Fig. 1).As a target species, three mammals of different taxa (European bison Bison bonasus, Eurasian red squirrel Sciurus vulgaris, and European hamster Cricetus cricetus), with different body size (0.2-800 kg), lifestyle (terrestrial, burrowing and arboreal species), social organization (living in herds or solitary) and different predicted type of hair glucocorticoids were selected.The selected species are protected in numerous European countries, and extensive research in the fields of ecology and conservation biology is being conducted for them.
We believe that the development of new methodological protocols for determining glucocorticoids in the hair of selected species will allow us to better understand the condition of individual populations and conduct their more effective conservation.Performing such studies on hair samples will allow us to assess long-term glucocorticoid levels and obtain important information about free-living animals with minimal impact on the study organisms.At the same time, we believe that the development of validated LC-MS/MS methods for various groups of wild animals will allow for the dissemination of such analyses in conservation biology.

Development of sample preparation protocol
Next, we adjusted various conditions of the sample preparation protocol to ensure the simultaneous extraction of these three important glucocorticoids in a single extraction step and their reliable quantification in hair samples from different animals.Preliminary studies showed that hair from European bison has the lowest content of studied glucocorticoids compared to other evaluated species (red squirrel, European hamster).Furthermore, it was possible to find a hair-blank sample (containing analytes at an undetectable level) in the amount required for the whole optimization process.The sample preparation step was optimized using hair extracts fortified with cortisol, cortisone, and corticosterone in appropriate concentrations alongside corresponding blank samples.
Prior extraction of glucocorticoids, hair samples from animals were washed twice with isopropanol.This allows for the removal of external contaminants (including cortisol deposited from sweat or sebum) and represents a standard step of hair preparation of many existing protocols 19,28 .However, methanol might be also used for this purpose 23 .To date, a double wash of hair shafts with isopropanol allows for the effective removal of external cortisol with minimal effect on its fraction within the internal hair matrix 19 .
The "gold-standard" method used for the extraction of different glucocorticoids (including cortisol, cortisone, and corticosterone) from the hair is overnight incubation of the sample with methanol 21,25,29,30 .Herein, this strategy was also employed.Conditions of the sample preparation protocol were adjusted for ~ 40 mg of the animal hair sample.Available LC-MS/MS protocols for steroid hormones usually require from 25 mg (horse) or 40 mg (red deer) to 1 g (polar bear) of hair for analysis 21,28,30 .
Two different strategies were evaluated to ensure adequate clean-up of the methanolic extract of animal hair before the determination of glucocorticoids: solid-phase extraction (SPE) and dispersive solid-phase extraction (dSPE).To this, the effectiveness of extract clean-up using four different commercially available dSPE sorbent mixtures and nine SPE cartridges providing different retention mechanisms was evaluated by determining extraction recoveries (R E ), apparent recoveries (ApR E ), and matrix effect (ME).Assessment of these parameters gives information about the potential loss of analytes during sample preparation (R E ), the influence of the co-extracted hair matrix component on the ionization of the analytes in ESI source (ME), and the overall method efficiency (ApR E ).ApR E s were calculated by dividing the glucocorticoid signal in a pre-extracted fortified blank hair sample and that in pure methanol.The obtained data are included in Fig. 2. Generally, noticeably better ApREs and less suppression of analyte ions were reported for SPE than dSPE purification.Regarding R E s and ApR E s, clean-up of hair extracts using SOLA HRP SPE cartridges and dSPE 4 (sorbent mixture containing MgSO 4 , primary secondary amine (PSA), graphitized carbon black (GCB), and C18) was unsuitable for the determination of glucocorticoids.Sample purification on STRATA-X (containing polymer-based SPE sorbent-a surface-modified styrene divinylebenzene) resulted in the best R E s for all studied analytes, and high ApR E s compared to other SPE cartridges.The type of SPE cartridges has been so far applied for the LC-MS/MS quantification of urinary cortisol, cortisone, and their tetrahydro-metabolites 20 .In agreement with other research data 20 , STRATA-X cartridges were found to be more adequate for analytes determination than HLB cartridges, which also contain modified styrene polymer as a packed bed.In the case of corticosterone, clean-up using Discovery MCAX, DCS-Ph, DCS-CN cartridges provided better ApR E s than Discovery C8, C18, and C18Lt once.Clean-up of hair extract using Discovery MCAX cartridge also resulted in better ApR E s for cortisol and cortisone compared to other tested SPE cartridges (except STRATA-X for both analytes, and Discovery C18Lt for cortisone) (Fig. 3).A series of experiments were then performed to determine the best conditions for extraction of glucocorticoids from hair followed by extract purification on STRATA-X SPE cartridges (Fig. 4).Analytes extraction from ground to fine powder hair shafts instead of cut once yielded ~ 40% and ~ 30-35% increase in R E s and ApR E s values, respectively.This observation is consistent with other studies 19 since grounding allows for the break up of the structure of the hair shaft and increases the surface area for extraction.Furthermore, the worst R E s for all analytes and slightly better ApR E s for cortisol and cortisone were reported when extraction was performed at 56 °C instead of room temperature (Fig. 4).Conducting extraction step at elevated temperature is recommended in some studies 24 .Also, a significant drop in R E s and ApR E s values was noted when the methanolic extract was directly loaded on the SPE cartridge without evaporation and sample reconstitution in water (aqueous extract).Notably, enhancement of glucocorticoids signals was noted for methanolic extracts (MEs > 100%), while for aqueous extracts signals suppression was observed (MEs < 80%).Based on the obtained results, we decided to pulverize the hair shafts before extraction (performed at room temperature) with methanol and to replace the sample solvent with water before SPE purification.
Additionally, we provided further improvements on the SPE step: (1) the elution solvent was changed from methanol/acetonitrile/3 mol/L acetic acid (40:40:20, v/v/v) to methanol, (2) the eluant was evaporated up to 200 µL instead of complete sample drying.These changes were made to shorten the evaporation step (⁓ 5 times) and avoid possible loss of analytes caused by sample overdrying or insufficient redissolving of the residue.Finally, we received slightly worse R E s (but still within acceptable ranges), and less ion suppression caused by hair matrix www.nature.com/scientificreports/for steroids, and better ApR E s for cortisone and corticosterone (Fig. 5).Furthermore, we reduced sample handling and sample preparation time before LC-MS/MS analysis.

Method validation
The developed method was validated using hair samples from three different animals (European bison, red squirrel, European hamster) to confirm its flexibility and practical utility.To improve both the accuracy and precision of the obtained results, data for each target compound was normalized on one internal standard-deuterated analog of cortisol (cortisol-D 4 ).The implementation of surrogate analytes (e.g.stable-isotope-labeled analogs, SILs) is a popular strategy for reliable quantification of glucocorticoids in different biological matrices including hair 18,24,25 .Undoubtedly, the main disadvantage of this approach is the high cost of SILs.Thus we verified if it is possible to obtain reliable results by applying one SIL for three glucocorticoids, to reduce cost per analysis and develop a protocol attractive for routine analysis.Depending on the evaluated animal species, the developed method allows for cortisol detection at 0.07-1.19ng/mL, cortisone at 0.05-0.46ng/mL, and corticosterone at 0.45-1.26ng/mL, which corresponds to 3.84 -29.77 pg/mg, 1.86-11.47pg/mg, and 11.31-31.51pg/mg in hair, respectively (Table 2).For all glucocorticoids, the best LODs in hair matrix were noted for the red squirrel and the worst for the European hamster.Furthermore, the method showed good linearity for all analytes in animal hair matrix regardless of species studied and a wide linear range of up to 5 ng/mg.Better LODs were reported in hair collected from the human scalp 24,25,31 or polar bear 28 .However, the obtained LOQs are sufficient to reliably monitor cortisol, cortisone, and corticosterone in the hair of the study group of free-living animals.
(Table 3).For European bison hair, the method showed good recoveries within the range of 83.37-104.97%for all glucocorticoids.For other species, in some cases, the exceedance of acceptable levels was observed (R E < 70%).At the sample preparation step, the most valuable loss was noted for corticosterone in red squirrel hair in quality control (QC) samples containing this analyte at low (LQC) and medium (MQC) concentration levels (R E s between 41.20 and 47.41%).In QC samples fortified with analytes at high concentration level (HQC), ⁓ 34-40% loss of corticosterone and cortisone were also reported for hair samples from European hamster and red squirrel, respectively (Table 3).Other researchers have reported better R E for corticosterone (~ 80%), and worse for cortisone (~ 68%) and cortisol (~ 55%) in human scalp hairs after pretreatment combining SPE and liquid-liquid extraction (LLE) with toluene 31 .Suppression of steroid hormone signals by hair matrix components has been widely reported 31 .In our study, co-extracted hair matrix components caused suppression of ionization (ME < 100%) for evaluated analytes, and the assessed ME can be classified as medium or strong (Table 3).
Only for cortisol determined in European hamster hair samples, MEs were negligible regardless of the studied concentration level.These results implied the complex composition of animal hair and confirmed the need for the additional purification of the methanolic hair extracts before LC-MS/MS analysis.Regarding ME data, ("cut hair") or further grounded using a standard procedure before extraction ("powdered hair").Next, analytes were extracted using 1 mL of methanol for 24 h at room temperature (see-saw rocker) or 56 °C (thermoblock).
The supernatants (800 µL) were loaded on SPE cartridge ("methanolic extract") or evaporated to dryness and reconstituted with 1 mL of water ("aqueous extract") before SPE purification.SPE was performed on STARTA-X cartridges and further treated as specified in Fig. 2. Experiments were carried out using blank samples and fortified with glucocorticoids at 0.05 µg/mL before extraction.www.nature.com/scientificreports/ the matrix-matched calibration was selected for quantitative purposes.Additionally, the analytes have shown satisfactory stability in final samples (S FT and S A values between 91.40 and 109.73%),which underlines practical usefulness of the developed method (Table 3).

Analysis of animal hair samples
Quantitative results of the concentration of cortisol, cortisone, and corticosterone in hair samples collected from European bison, European hamster, and Eurasian red squirrel are summarized in Table 4 and supported by the examples of LC-MS/MS data in Supplementary Fig. 1S.The positive values obtained for the tested glucocorticoids ranged from 5.7 to 215.2 pg/mg, however, ND (i.e., "not detected") cases were also obtained.In most cases, a "not detected" result is expected and is not related to the quality of the method.It is known that in certain animal species, the metabolic pathway of cortisol (in European bison and red squirrel) or corticosterone (in European hamster) is dominant.In this context, the lack of e.g., corticosterone in the hair of European bison is a correct and expected result.However, "not detected" results for the cortisol-cortisone pair of hormones in animals where cortisol is the primary glucocorticoid should be interpreted differently.It is known that cortisol can be reversibly transformed to its inactive metabolite-cortisone and the detection of both steroids is likely.For example, in the red squirrel, both hormones were found in significant concentrations in every sample tested.In the European bison hair, however, either one or the other glucocorticoid (more often cortisone) was detected (Table 4; Fig. 6).If only one of the hormones (cortisol or cortisone) were utilized for the research, entirely different animals would be identified as potentially the most susceptible to the analyzed stressors.The obtained results for European bison show how important it is for the correct assessment of stress reactions of protected animal species to simultaneously determine several glucocorticoids from a given metabolic pathway.Our method clearly shows that measuring only cortisol could lead to erroneous conclusions about which individual has higher levels of stress hormones.Such a result may be an important clue in studies that use only enzyme immunoassays and determine only a single compound.As our method shows, at least for some species the set of analyzed makers should be extended.

Optimization of LC-MS/MS parameters
In initial conditions, the scan spectra of target analytes (including internal standard) were acquired in ESI + to select the precursor ions.ESI in positive polarity allows for the detection of low levels of glucocorticoids using LC-MS/MS in different biological matrices 18,21,24 .The ions corresponding to [M + H] + were selected to develop the quantitative method.Next, the fragmentor voltage (60-160 V) and collision energy (0-40 eV) were optimized to provide appropriate conditions for the formation of precursor ions.The obtained data were inspected to select characteristic transitions and selectively detect analytes in a multiple reaction monitoring (MRM) mode.Next, key parameters of the mass spectrometer (MS) were optimized to achieve good sensitivity of measurement: nebulizer pressure, drying gas temperature, drying gas flow, sheath gas temperature, sheath gas flow, capillary voltage, and nozzle voltage within the range of 30-45 psi, 250-350 °C, 5-9 L/min, 250-400 °C, 5-10 L/min, 3000-5000 V, and 0-900 V, respectively.

Hair sample preparation
Hair strands were washed twice by hand shaking in 5 mL of isopropanol for 3 min to remove the extrinsic contamination from their surface.The washing solutions were decanted and removed, and hair samples were dried overnight at room temperature.After this step, the hair was carefully cut with surgical scissors into ⁓ 5 mm pieces.Then, ~ 40 mg of hair sample was accurately weighed using the OHAUS Adventrer (Parsippany, New Jersey, USA), transferred into a 2 mL plastic tube, and milled at 3000 rpm with five steel beads (3 mm) for 9 min (in three cycles of 3 min each).The powdered hair sample was incubated in 1 mL of methanol containing 10 µg/mL of cortisol-D 4 at room temperature using a Stuart SSL4 see-saw rocker (Vernon Hills, Illinois, USA) at 70 rpm speed range.After 24 h, the samples were centrifuged (13 000 rpm, 15 min) using a 5415R centrifuge (Eppendorf, Germany), and supernatants (800 µL) were transferred to a new centrifugal tube.Then, aliquots www.nature.com/scientificreports/were evaporated to dryness using an EZ-2 Elite Personal Evaporator (Genevac Ltd., Ipswich, UK), reconstituted in 1 mL of ultrapure water, and sonicated for 5 min in an ultrasonic bath (Polsonic, Warsaw, Poland).Next, the obtained samples were purified by SPE using a Gilson GX-271 ASPEC (Middleton, WI, USA) automatic extraction system operated by Gilson Ethernet Utility 1.8.6.1 and Trilution LH 2.0 software.The Strata-X SPE cartridges were conditioned with 1 mL of methanol, and equilibrated with 1 mL of water, respectively.Then 1 mL of the sample was loaded onto the cartridge, which was further washed with 1 mL of water.Elution was carried out with 1 mL of methanol.The eluate was dried up to 200 µL.The clear supernatant was transferred into a glass chromatographic insert vial and analyzed in triplicate by the LC-MS/MS method.The injection volume was 10 µL.

LC-MS/MS analysis of hair extracts
Measurements were carried out on a 1290 infinity UHPLC system consisting of a degasser, binary pump, autosampler, and column thermostat (Agilent Technologies, Santa-Clara, CA, USA  1.
An amount of individual glucocorticoid in authentic animal hair samples was determined from the corresponding matrix-matched calibration curve.Next, the values in µg/mL were converted to µg/g considering the mass of hair strands subjected to milling and extraction.The data were expressed as the mean value from three injections of the sample plus the standard deviation.

Method validation
For validation experiments, representative hair pools originating from different individuals showing undetectable (from European hamster and European bison) or low basal levels (for red squirrel) of target glucocorticoids.A series of quality control (QC) samples were prepared to evaluate: linearity, the limit of detection (LOD), the limit of quantification (LOQ), accuracy, precision, extraction recovery (R E ), autosampler stability (S A ), the freeze-thaw stability (S FT ), and matrix effect (ME).Blank (unspiked samples) were also analyzed, and the peak area corresponding to the endogenous concentration of the analyte in the hair pool was subtracted from the peak area detected in the QC sample.Acceptance criteria were taken from guidance published by FDA 32 .
For linearity evaluation, nine matrix-matched calibration standards with increasing concentrations of glucocorticoids in the range of 0.3 ng/mL to 0.2 µg/mL were analyzed.Each sample was fortified with a constant amount of internal standard.Next, the regression plots were built for each analyte using the response factorsthe ratio of the analyte peak area over cortisol-D 4 peak area.Linearity was acceptable with the coefficient of determination (R 2 ) of more than 0.99.For each animal species, LODs and LOQs were calculated independently using the formula: LOD or LOQ = a × SD/b, where a is a factor of 3.3 or 10 for LOD and LOQ, respectively, SD is a standard deviation of the intercept and b is a slope of the regression line.
Accuracy, precision, stability, R E and ME were studied for hair matrix using QC samples fortified with analytes at three different concentration levels: low (LQC = ~ 3 × LOQ), medium (MQC), and high (HQC).Interday and intraday variations of accuracy (trueness) and precision (repeatability) were calculated using data obtained from 3 or 7 individually prepared QC samples analyzed within one day or on three different days, respectively.The precision was expressed as the relative standard deviation (RSD), while accuracy was calculated using ratios of analyte concentrations determined in the samples from the corresponding calibration curve (A) to concentration added (B) (accuracy = A/B × 100%).Acceptable levels of accuracy and precision were within 15%.For R E , the ratios of peak areas in pre-extraction (C) and post-extraction (D) spiked pooled hair matrix were compared (R E = (C/D) × 100%; acceptable range: 70-120%).For ME calculation, the ratios of peak areas in post-extraction spiked hair sample (D) to corresponding spiked neat solvent (E) were compared (ME = (D/E) × 100%).ME values within 0-10%, 10-20%, 20-50%, and < 50% were classified as negligible, soft, medium, and strong, respectively.The stability was estimated by comparing signals of glucocorticoids in freshly prepared hair extracts to samples stored for one day at 25 °C in the autosampler tray (for S A ) or samples subjected to three freeze/thaw cycles from 25 to − 20 °C (for S FT ).Acceptable levels of stability were within ± 10%.Each QC sample was injected three times.

Analysis of hair samples from free-living animals
To demonstrate the applicability of the developed analytical method, the determinations of glucocorticoids in hair samples from free-living animals of selected species were carried out.We chose three mammal species with a high conservation status for our study.Protected species and their welfare are frequent subjects of research, making the development of dedicated analytical methods appealing to a broader audience.We selected animals from various taxonomic groups: European bison, European hamster, and Eurasian red squirrel.These animals differed in terms of body size, lifestyle (including terrestrial, burrowing, and arboreal species), social organization (either living in herds or solitary), and the expected type of hair glucocorticoids.
European bison is a symbolic species for nature conservation in Europe.It is the largest living mammal in Europe.This species became extinct in the wild at the beginning of the twentieth century, but thanks to many conservation efforts, its population has been restored and currently there are about 10,000 European bison in www.nature.com/scientificreports/Europe, including nearly 80% in the wild.This species occurs mainly in forests, but it also willingly uses open areas, which provide it with an attractive food base [33][34][35] .Due to various threats [36][37][38] , European bison populations are constantly monitored and the measurement of glucocorticoids can be a good tool that can be routinely used in monitoring.The development of an analytical method for measuring corticoids in hair for this species is therefore very justified from a practical point of view.European hamster is a burrowing rodent species with high conservation status in the European Union.It is strictly protected in many European countries, listed in EU legislation (Habitats Directive) and international conventions (Appendix II of the Bern Convention).This species is totally dependent on agriculture because farmlands and human settlements occupy its original fertile steppe habitats 39 .Nowadays, many research and conservation actions are being taken for this mammal [40][41][42][43][44][45] .For this reason, the measurement of glucocorticoids can be a useful tool that can be routinely used in the monitoring of endangered or reintroduced hamster populations.
The Eurasian red squirrel is also a protected species in Poland.Due to the occupation of urbanized areas, this species is exposed to various anthropogenic stressors 46 .At the same time, in some European countries, this species is under pressure from the geographically alien (originating from North America) eastern grey squirrel Sciurus carolinensis, which is also a carrier of the highly pathogenic to red squirrels parapoxvirus 47 .It is also known from the literature that squirrels are characterized by quite high levels of steroids [48][49][50] and thus can be a very useful species to test analytical methods as a reference point for other species.

Conclusions
An accurate LC-MS/MS method for the determination of cortisol, cortisone, and corticosterone in animal hair was developed and validated.The method flexibility was evaluated on hair samples of three different species of mammals with high conservation status and different ecological features.The method comprises a simple extraction with methanol and further sample clean-up using solid-phase extraction (SPE).A special emphasis was put on the purification of hair extracts to reduce the suppression of the glucocorticoid's signals caused by the matrix background.Numerous commercially available SPE cartridges and dispersive SPE kits were selected for hair extract purification.The method performance was evaluated on representative hair pools delivered from selected mammal species by determining key validation parameters.
With increasing human pressure and climate change, there is an urgent need to develop validated tools for reliable assessment and monitoring of the environmental stressors for various free-living animals.This allows for a better understanding of contemporary natural processes and taking effective remedial action.The method presented in this work demonstrates that keratinized tissues of animals can be successfully utilized for research, and measuring multiple glucocorticoids simultaneously provides more comprehensive information than solely concentrating on one of them.However, the utility of specific glucocorticoids as markers of the stress response in keratinized tissues should be substantiated by additional experimental studies on specific animal targets.During these studies, possible pathways of glucocorticoid deposition in hair, as well as species differences in the hair growth and shedding cycle, should be considered.Attention should also be given to the potential specificity of glucocorticoid deposition in the hair of hibernating species, where samples collected in spring may reflect the free corticosterone concentration from the hair growth period before hibernation.

Figure 1 .
Figure 1.Scheme summarizing the crucial aims of the work.

Figure 4 .
Figure 4. Effect of different strategies of analytes extraction and extract pretreatment on extraction recovery (A), matrix effect (B), and apparent recovery (C).A representative sample of European bison hair was cut into small segments with scissors, and transferred into tubes (40 mg), and directly subjected to extraction ("cut hair") or further grounded using a standard procedure before extraction ("powdered hair").Next, analytes were extracted using 1 mL of methanol for 24 h at room temperature (see-saw rocker) or 56 °C (thermoblock).The supernatants (800 µL) were loaded on SPE cartridge ("methanolic extract") or evaporated to dryness and reconstituted with 1 mL of water ("aqueous extract") before SPE purification.SPE was performed on STARTA-X cartridges and further treated as specified in Fig.2.Experiments were carried out using blank samples and fortified with glucocorticoids at 0.05 µg/mL before extraction.

Figure 5 .
Figure 5.Comparison of extraction recovery (A), matrix effect (B), and apparent recovery (C) obtained during testing of different protocols for elution of glucocorticoids from SPE cartridges and sample concentration before LC-MS/MS analysis.Experiments were conducted using 40 mg of grounded European bison hair sample extracted with 1 mL of methanol (containing 0.05 µg/mL of cortisol, cortisone, and corticosterone) for 24 h at room temperature.Supernatants (800 µL) were evaporated to dryness and redissolved in 1 mL of methanol before being loaded on STRATA-X SPE cartridges.During SPE clean-up, analytes were eluted with 1 mL of methanol/acetonitrile/3 mol/L acetic acid (40:40:20, v/v/v) (procedure 1) or methanol (procedure 2).Next, 1 mL of the sample was evaporated to dryness and redissolved in 200 µL of methanol (ultrasonic-assisted step, 2 min)-procedure 1, or its volume was reduced up to 200 µL by evaporation-procedure 2. https://doi.org/10.1038/s41598-023-49940-2

Table 2 .
Comparison of limits of detection (LODs) and limits of quantification (LOQs), and key parameters of matrix-matched calibration curves estimated for glucocorticoids in European bison, Eurasian red squirrel, and European hamster hair matrices.LR linear range of calibration curve, REq regression equation, R 2 coefficient of determination.*x in µg/mL; z in ng/mg.

Table 4 .
The The hair samples used for this study came from a wild animal tissue bank existing at the following universities in Poland (central Europe): Warsaw University of Life Sciences (consent of Regional Directorate for Nature Protection in Warsaw for red squirrel number WPN-I.6401.208.2018.PF, and for European bison number WPN-I.6401.90.2014.EB.1) and the John Paul II Catholic University of Lublin (consent of Regional Directorate for Nature Protection in Lublin number: WPN.6401.3.1.2023.AO).No animal was intentionally captured or killed for this work.
concentration of cortisol, cortisone, and corticosterone obtained during LC-MS/MS analysis of experimental hair samples from free-living European bison, European hamster, and Eurasian red squirrel.ND "not detected"; D detected, not quantified.Cortisol Cortisone Figure 6.Concentration of cortisol and cortisone in hair samples from 10 European bison individuals.Corticosteroids were not consistently present in all samples, and the concentration of cortisol did not correlate with cortisone.Vol:.(1234567890)Scientific Reports | (2023) 13:23089 | https://doi.org/10.1038/s41598-023-49940-2www.nature.com/scientificreports/Hair samples characteristics

Table 5 .
Characteristics of the samples used for analysis include information on the animal species from which the hair samples were collected, the animal's sex, age or body weight, its approximate location and the date of sampling.

mass The location of the sample's origin Date of sample collection
).The UHPLC instrument set-up was connected with an Agilent 6460 triple quadrupole mass spectrometer (QQQ) equipped with an Agilent Jet Stream electrospray ion (ESI) source.Data were acquired and analyzed with Agilent MassHunter Acquisition software v.B.08 and Agilent MassHunter Quantitative Analysis software v.B.07, respectively.Chromatographic separation was achieved at 40 °C on the Zorbax Eclipse Plus C18 RRHT column (2.1 mm × 100 mm × 1.8 µm) protected by the Zorbax Eclipse Plus C18 guard column (2.1 mm × 12.5 mm × 1.8 µm) both from Agilent Technologies.Gradient elution was performed using a mobile phase consisting of 0.1% (v/v) formic acid in water (solvent A) and methanol (solvent B) at the flow rate of 0.35 mL/min and the following program: 0-4 min 20-60% solvent B; 4-5 min 60% solvent B; 5-7 min 60-70% solvent B (post run: 2 min).The MS conditions were optimized using Source and iFunnel Optimizer software v.B.08.00 (Agilent Technologies).Data were acquired in dynamic multiple reaction monitoring (dMRM) mode in positive polarity using settings detailed in Table